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The  contents  of  this  report  reflect  the  views  of  the  Jet 
Propulsion  Laboratory  which  is  responsible  for  the  facts 
and  accuracy  of  data  presented.  This  report  does  not  con- 
stitute a standard,  specification  or  regulation. 
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A major  mission  of  the  U.S.  Coast  Guard  is  the  task  of  providing  and  maintaining  Mari- 
time Aids  to  Navigation.  These  aids  are  located  on  and  near  the  coastline  and  inland 
waters  of  the  United  States  and  its  possessions.  A computer  program.  Design  Synthesis 
and  Performance  Analysis  (DSPA),  has  been  developed  by  the  Jet  Propulsion  Laboratory 
to  demonstrate  the  feasibility  of  low-coat  solar  array/battery  power  systems  for  use 
on  flashing  lamp  buoys.  ^ 

These  buoys  are  ejcp&sSdi  to  a wide  variety  of  environments.  The  range  of  these  environ- 
ments are  typified  by  weather  conditions  in  locations  in  Alaska  and  the  Florida  coastal 
waters..*?’To  provide  detailed,  realistic  temperature,  wind,  and  solar  insolation  data 
for  analysis  of  the  flashing  lamp  buoy  power  systems ^the> Jet  Propulsion , Laboratory 
developed* two  DSPA  support  computer  program  sets:  MERGE  and  STAT.'  A' general 
description  of  these  two  packages  is  presented  in  this  program  summary  report. 

The  MERGE  program  set  will  enable  the  Coast  Guard  to  combine  temperature  and  wind  velo- 
city data  (NOAA  TDF-lU  tapes)  with  solar  insolation  data  (NOAA  DECK-280  tapes)  onto  a 
single  sequential  (^MERGE'^^ile  containing  up  to  12  years  of  hourly  observations.  This 
MERGE  file  can  then  be  used  as  direct  input  to  the  DSPA  program. 

-O  - 

The  STAT  program  set  will  enable  th»-Caaat  Guard  "t'o  pei’fero^a  statistical  analysis^of 
the  MERGE  data  and  produce  high  or  low  or  mean  profiles  of  the  data  and/or  do  a worst 
case  analysis.  The  STAT  output  file  consists  of  a one-year  set  of  hourly  statistical 
weather  data  which  can  be  used  as  input  to  the  DSPA  program. 
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INTRODUCTION  AND  SUMMARY 


A major  mission  of  the  U.S.  Coast  Guard  is  the  task  of  providing  and 
maintaining  liaritime  Aids  to  Navigation.  These  aids  are  located  on  and 
near  the  coastline  and  inland  waters  of  the  United  States  and  its  pos- 
sessions. A computer  program.  Design  Synthesis  and  Performance  Analysis 
(DSPA),  has  been  developed  by  the  Jet  Propulsion  Laboratory  to  demonstrate 
the  feasibility  of  low-cost  solar  array /battery  power  systems  for  use  on 
flashing  lamp  buoys. 

These  buoys  are  exposed  to  a wide  variety  of  environments.  The  range  of 
these  environments  are  typified  by  weather  conditions  in  locations  in 
Alaska  and  the  Florida  Coastal  waters.  To  provide  detailed,  realistic 
temperature,  wind,  and  solar  insolation  data  for  analysis  of  the  flashing 
lamp  buoy  power  systems,  the  Jet  Propulsion  Laboratory  developed  two  DSPA 
support  computer  program  sets:  IlERGE  and  STAT.  A general  description  of 
these  two  packages  is  presented  in  this  program  summary  report.  Details 
on  the  content  and  usage  of  these  packages  are  presented  in  the  following 
reports : 

Program  Documentation  (DS/PA  Computer  Program): 

(JPL  Report  No.  50U0-27) 

Volume  I - Software  Requirements  Document. 

Volume  II  - User's  Manual. 

Volume  III  - Programmer's  Manual. 

The  MERGE  program  set  will  enable  the  Coast  Guard  to  combine  temperature 
and  wind  velocity  data  (NOAA  TDF-lU  tapes)  with  solar  insolation  data 
(NOAA  DECK-280  tapes)  onto  a single  sequential  "ivERGE"  file  containing  up 
to  12  years  of  hourly  observations.  This  MERGE  file  can  then  be  used  as 
direct  input  to  the  DSPA  program. 

The  STAT  program  set  will  enable  the  Coast  Guard  to  perform  a statistical 
analysis  of  the  MERGE  data  and  produce  high  or  low  or  mean  profiles  of 
the  data  and/or  do  a worst  case  analysis.  The  STAT  output  file  consists 
of  a one-year  set  of  hourly  statistical  weather  data  which  can  be  used 
as  input  to  the  DSPA  program. 


i 


The  final  phase  of  the  development  of  the  DSPA  support  programs  was  the  ^ 
production  of  MERGE  and  intermediate  STATS  (MERGE  averages  file  for  sub-  ; 
sequent  profiling)  files  for  each  of  the  15  Coast  Guard  selected  loca-  : 
tions.  The  NOAA  TDF-lU  and  DECK-280  data  used  for  generating  these  ; 
files  are  shown  in  Table  1-1  below.  : 
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TABLE  1-1. -INDEX  TO  NOAA  SOLAR/WEATHER  TAPES 


Reel 

No. 

Station 

Number 

Station  Location 

Time  Period  Covered 

X788 

12839 

Miami,  Florida 

Jan.  1955  to  Dec.  1964 

X789 

12919 

Brownsville,  Texas 

Jan.  1953  to  Dec.  1956 
Jan.  1959  to  Dec.  1964 

X790* 

13745 

Hatteras,  North  Carolina 

Jan.  1955  to  Feb.  1957 

(93729) 

(Cape  Hattsras,  North 

Carolina 

Mar.  1957  to  Dec.  1965 

X791 

14607 

Caribou,  Maine 

Jan.  1955  to  Dec.  1964 

X792 

14732 

La  Guardia  Field,  New  York 

Jan.  1953  to  Dec.  1958 
Jan.  1951  to  Dec.  1954 

X793 

14739 

Boston,  Massachusetts 

Jan.  1955  to  Dec.  1964. 

X794 

14847 

Saulte  Ste.  Marie, 

Michigan 

Jan.  1952  to  Dec.  1961 

X795 

23174 

Los  Angeles,  California 

Jan.  1955  to  Dec.  1964 

X796 

24233 

Seattle,  Washington 

Jan.  1955  to  Dec.  1964 

X797 

26615 

Bscrel,  Alaska 

Jen.  1952  to  C-3C.  T-;: 

X799 

13743 

Washington,  D.C. 

Jan.  1955  to  Dec.  1964 

X800 

13983 

Columbia,  Missouri 

Jan.  1955  to  Dec.  1964 

X801 

23154 

Ely,  Nevada 

Jan.  1955  to  Dec.  1964 

X802 

24225 

Medford,  Oregon 

Jan.  1955  to  Dec.  1554 

X803 

93193 

Fresno,  California 

Jan.  1955  to  Dec.  1964 

12839 

Miami,  Florida 

Jan.  1955  to  Dec.  1964 

12919 

Brownsville,  le'^as 

Jan.  1953  to  De:.  ''555 
Jan.  1 959  to  Dec.  1 555 

13745 

Hatteras,  North  Carolina 

Jan.  1955  to  Feb.  1957 

14607 

14753 

Caribou,  Maine 

Blue  Hill,  Massachusetts 

Jan.  1955  to  Dec.  1964 
Jan.  1955 

X305** 

14753 

Blue  Hill,  Massocnusests 
(cont) 

to  Dec.  1954 

14847 

Saulte  Ste.  Marie, 

Michigan 

July  1952  to  Aug.  1953 

23174 

Los  Angeles,  California 

Jan.  1962  to  Dec.  1966 

24233 

Seattle,  Washington 

Jan.  1955  to  Dec.  1564 

26615 

Bethel,  Alaska 

July  1952  to  Dct.  1552 
Dec.  1956  to  Apr.  1957 

X806** 

93729 

Capa  Hatteras,  North 

Carolina 

Mar.  1957  to  Dec.  1364 

94706 

New  York,  New  York 

Jan.  1953  to  Dec.  1958 
Jan.  1963  to  Dec.  ).'f4 

X807 

13933 

Columbia,  Missouri 

Jan.  1955  to  Dec.  1964 

23154 

Ely,  Nevada 

Jan.  1955  to  Dec.  1964 

XS08 

24225 

Medford,  Oregon 

Jan.  1955  to  Dec.  1564 

93193 

Fresno,  California 

Jan.  1955  to  Dec.  1964 

X809 

93722 

Silver  Hill,  Maryland*** 

Jan.  1955  to  Dec.  1560 

93734 

Sterling,  Virginia*** 

Jan.  1961  to  Dec.  1964 

* The  Hattoras  and  Cape  Hatteras TDF14  d’i.a  '.'.d  to  te  ccr pined  into  one  file  to 

Facilitate  the  building  of  the  sleletcn  Hir/jii  iile 

**  No  End-of-File  mark  exists  between  jtatlons  - this  problem  Is  handled  by  the 

DECK280  program. 

***  Solar  insolation  data  from  these  two  files  arc  for  use  with  the  Washington,  O.C. 
TDF14  data. 
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2.  MERGE  COMPUTER  PROGRAM  SET  DESCRIPTION 


The  features  of  the  MERGE  computer  program  set,  consisting  of  the  TDFll*, 
DECK280,  and  LISTMERGE  routines,  are  discussed  in  this  section.  The 
elements  of  this  discussion  include  specific  program  requirements, 
computer  program  rationale,  TDFlU  program,  DECK280  program,  and 
LISTMERGE  program.  These  elements  are  presented  below. 

2.1  Specific  Program  Requirements 

The  major  requirement  of  the  MERGE  computer  program  set  is  to  combine  the 
solar  insolation  data  from  the  NOAA  DECK-280  tapes  with  the  temperature 
and  wind  velocity  data  from  the  TDF-ll*  tapes.  The  resultant  output,  a 
MERGE  file,  is  to  be  a computer  data  file  consisting  of  no  more  than 
12  years  of  information.  The  MERGE  file  information  is  to  be  recorded  on 
an  hourly  basis  and  stored  on  one-day  records,  all  days  being  stored  in 
Julian  date  sequence.  Temperature  data  is  in  degrees  Fahrenheit j wind 
velocity  is  in  knots,  and  solar  insolation  is  in  watts/square  meter. 


2.2  Computer  Program  Rationale 

The  basis  for  selecting  the  proper  computer  program  for  meeting  the  Coast 
Guard  requirements  was  determined  by  examination  of  the  necessary  DSPA 
weather  data  input.  Generally  the  DSPA  program  requires  that  temperature 
and  solar  insolation  values  be  known  (wind  velocity  was  added  as  a 
possible  future  i-equirement ) . Accuracy  requirements  further  dictated 
that  data  observations  should  be  hourly.  Use  of  the  NOAA  weather  tapes 
provided  the  needed  data,  though  not  in  a very  accessible  form. 

Investigation  of  the  NOAA  tapes  revealed  that  the  TDF-lU  data  was  edited 
and  sequentially  stored,  but  the  DECK-280  data  was  unedited  and  unsorted. 
Hence,  it  was  determined  that  the  MERGE  process  should  be  performed  in 
two  separate  steps: 

1)  Generation  of  a skeletal  MERGE  file  from  the  TDF-lU  tape  data. 

2)  Addition  of  the  DECK-280  data  to  the  MERGE  file. 
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An  overview  flow  diagram  of  the  MERGE  program  ia  ohown  in  Figure  a-i  below. 
A sample  of  the  MERGE  output  file  contents  is  provided  in  Figure  2-2. 
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FIGURE  2-2.  SAMPLE  MERGE  OUTPUT  (Sheet  2 of  S) 


enter  start  o*te  ittOdd)  and  NOt  Or  Records  to  disPlat 


FIGURE  2-2.  SAMPLE  MERGE  OUTPUT  (Sheet  3 of  3) 
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I . 3 TDFlt4  Program 

The  TDFlU  program,  shown  in  Figure  2-3,  controls  the  building  of  the 
skeletal  MEHOE  file  from  the  IJOAA  TDF-lt  data  tape.  The  user  requests 
the  creation  of  a MERGE  file  to  span  a particular  periou  of  years  for  a 
selected  TDF-lt  location.  The  TDFlU  program  then  extracts  the  d 
tempei'ature , and  wind  velocity  information  from  the  input  tape.  ds 

a one-day  record  consisting  of  2b  hourly  observations  of  temperatur-'  and 
wind  velocity  and  space  for  solar  insolation,  and  sequentially  writer 
the  day's  data  to  a MERGE  file. 


FIGURE  2-3.  TDFlU  COMPUTER  PROGRAM 


2.b  DECK280  Program 


Ad  iition  of  solar  insolation  data  to  a MERGE  file  created  by  TDFlb  is 
periormed  by  the  DECK230  program  (Figure  2-U).  The  user  requests  that, 
for  a particular  MERGE  location,  NOAA  DECK-280  tape  data  from  a specified 
location  be  inserted  into  the  file.  The  DECK280  program  extracts  the 
solar  radiation  data  (in  Langleys)  from  the  NOAA  tape,  converts  the  data 


FIGURE  2-k.  DECK280  COMPUTER  PROGRAM 
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to  watts/square  meter,  and  adds  the  data  to  the  appropriate  day  and  hour 
position  in  the  MERGE  file. 


2,5  LISTMERGE  Program 


The  average  MERGE  file  consists  of  from  10  to  12  years  of  hourly  tempera- 
ture, wind  velocity,  and  solar  insolation  data.  The  LISTMERGE  program 
permits  the  user  to  randomly  view  any  number  of  sequential  days  beginning 
at  any  date  contained  within  the  file.  (See  Figure  2-5). 


FIGURE  2-5.  LISTMERGE  COMPUTER  PROGRAM 
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\ 3.  STAT  COMPUTER  PROGRAM  SET  DESCRIPTIOH  j 

^ The  features  of  the  STAT  computer  program  set,  consisting  of  the  STATS 

‘ and  PROFILE  routines,  are  discussed  in  this  section.  The  elements  of 

i this  discussion  include  specific  program  requirements,  computer  program 

[ rationale,  STATS  program,  and  PROFILE  program.  These  elements  are  pre- 

f 

! sented  below. 

3.1  Specific  Program  Requirements 

The  major  requirement  of  the  STAT  computer  program  set  is  to  perform 
a statistical  analysis  of  the  MERGE  weather  data,  producing  high  or  low 
or  mean  or  worst  case  environmental  profiles.  The  user  must  be  able  to 

control  the  selection  of  alternative  profile  types  on  a monthly  basis. 

The  resultant  output  is  to  be  a computer  data  file  consisting  of  one  year 
of  statistical  weather  information  stored  in  the  same  format  as  the  tlERGE 
file  data.  j 

3 . 2 Computer  Program  Rationale  • 

The  basis  for  selecting  the  proper  computer  program  for  meeting  the  Coast  j 

Guard  requirements  was  determined  by  examination  of  the  process  involved  ■ 

in  profiling  data.  To  produce  a profile,  the  MERGE  data  must  first  be 
averaged  over  all  years.  This  "means"  data  represents  a "smoothed  over" 
year  with  an  extremely  small  likelihood  of  occurrence.  The  user,  there- 
fore, needs  to  be  able  to  adjust  the  information  by  applying  some  scaling 
factor  to  the  averaged  data  and,  if  desired,  performing  a worst  case 
analysis  of  the  data.  The  Coast  Guard  required  that  the  user  be  permitted 
to  treat  each  month  separately,  yielding  different  scale  factors  and/or 
worst  case  periods. 

Investigations  indicated  that  a process  Involving  reading  and  averaging 
and  profiling  of  the  10  to  12  years  of  MERGE  data,  all  in  one  step  and  for 
each  time  a new  statistical  profile  was  required,  would  be  prohibitively 
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expensive.  Hence,  it  was  determined  that  the  statistical  analysis  process 

should  be  performed  in  two  separate  steps: 

1)  Generation  of  a one-year  averages  file,  directly  from  the  MERGE  file, 
for  permanent  storage  on  r.agnetic  tape  (i.e.,  executed  one  time  only 
for  each  location). 

2)  Profiling  of  the  averages  data,  as  required,  for  input  to  th  "PA 
program  (approximately  4 times  less  expensive  than  step  1). 


An  overview  flow  diagram  of  the  ETAT  program  is  shown  in  Figure  3-1 
below.  Samples  of  the  STATS  and  PROFILE  output  file->  are  provided  in 
Figures  3-2  and  3-3  respectively. 


FIGURE  3-1  STAT  COMPUTER  PROGRAil  SET  OVERVIEW/ 


3. 3 STATS  Program 

The  STA.TS  program  uses  the  10  to  12  years  of  MERGE  file  weather  data  to 
produce  a one-year  statistical  file.  The  procedure  involves  the  averag- 
ing of  data  for  a given  hour  of  a given  day  for  each  of  the  years  con- 
tained in  the  MERGE  file.  The  statistical  data  is  then  written  onto  a 
STATS  file  which  is  used  as  input  to  the  PROFILE  program  (Section  3.4 
below).  Specifically,  the  STATS  program  computes  and  outputs  (both 
to  the  STATS  and  to  the  printer)  the  following  statistics: 

la)  Average  temperature  for  each  hour  of  one  year. 

lb)  Average  wind  velocity  for  each  hour  of  one  year. 

lc)  Average  solar  insolation  for  each  hour  of  one  year. 

2a)  Average  wind  velocity  for  each  day  of  each  data  year. 

2b)  Average  solar  insolation  for  each  day  of  each  data  year. 
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FIGURE  3-2B.  SAMPLE  STATS  OUTPUT-MONTHLY  DATA  (Sheet  2 of  2) 


FIGURE  3-2C.  SAMPLE  STATS  OUTPUT -YEARLY  DATA 
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FIGURE  3-3.  SAMPLE  PROFILE  PROGRAM  OUTPUT  (Sheet  3 of  6) 
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FIGURE  3-3.  SAMPLE  PROFILE  PROGRAM  OUTPUT  (Sheet  1*  ol  6) 
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3a) 

Average  ; 

3b) 

Average  ; 

ta) 

Average  " 

4b ) 

Average  \ 

4c ) 

Average  : 

5a) 

Average  i 

5b) 

Average  \ 

5c) 

Average  i 

6a) 

Standard 

6b ) 

Standard 

6c ) 

Standard 

7a) 

Minimum  i 

7b) 

Maximum  ■ 

8a) 

Mean  and 

8b) 

Mean  and 

For  statistics  la,^ lb,  and  Ic  the  STATS  program  averages  the  MERGE  values 
stored  foi?  a given  hour,  day,  and  month  over  all  data  years  contained  in 
the  file.  The  total  number  of  values  computed  and  stored  for  la,  lb,  or 
Ic  is  thus  2k  X 365  = 8t60,  giving  one  complete  year  of  average  hourly 
temperature,  wind  velocity,  and  solar  insolation  values.  This  constitutes 
the  bulk  of  the  STATS  file  content. 


These  hourly  averages  can  be  used  to  give  a reasonable  weather  profile 
for  any  hour  of  a sample  year.  However,  using  mean  values  for  an  entire 
year  produces  an  uncharacteristically  mild  profile  (a  mean  year)  with  an 
extremely  small  likelihood  of  occurrence.  Statistics  2a  through  8b,  as 
described  below,  are  therefore  gathered  to  allow  the  user  to  scale  the 
hourly  average  values  up  or  down  (Section  3.^  discusses  this  scaling). 

For  each  data  year,  the  STATS  program  averages  the  24  hourly  values  for  a 
given  day  to  produce  statistics  2a  and  2b.  The  total  number  of  values 
computed  and  stored,  therefore,  is  365  x for  each  statistic  (where 

is  the  number  of  years  of  data  available  for  location  1) . 
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The  STATS  program  then  averages  the  values  of  2a  or  2b  to  get 
statistics  3a  and  3b,  respectively,  for  each  day  of  the  year.  Thus, 
the  total  number  of  values  stored  for  each  of  these  statistics  is  365. 

For  statistics  Ua,  Ub,  and  Ac,  averages  of  all  hourly  temperature,  wind 
velocity,  and  solar  insolation  values,  respectively,  are  computed 
corresponding  to  a given  month  of  each  data  year.  The  total  n imber  of 
stored  values  for  each  of  the  Aa,  Ab,  or  Ac  statistics  is,  therefore, 

12  X 

The  program  averages  the  values  of  Aa,  Ab,  or  Ac  to  generate  statistics 
5a,  5b,  or  5c,  respectively,  and  computes  the  standard  deviation  (6a,  6b, 
or  6c)  of  the  Y^^  values  corresponding  to  each  month  of  the  year.  Hence, 
the  total  number  of  values  stored  for  each  of  these  six  statistics  is  12. 

Statistics  7a  and  7b  are  simply  the  minimum  and  maximum  temperatures 
observed  at  a location  during  a given  data  year.  Thus,  there  are  Y^ 
values  of  each  of  these  statistics. 

Finally,  the  program  computes  the  means  and  standard  deviations  of  the 
Y^  values  of  7a  or  7b  to  yield  statistics  8a  or  8b,  respectively.  Each 
of  the  8a  and  8b  statistics  has,  therefore,  only  two  stored  values. 

The  mathematical  formulas  for  the  above  statistics  are  provided  in 
Appendix  A.  A flow  diagram  of  the  STATS  program  is  shown  in  Figure  3-A. 


The  PROFILE  program  uses  the  one  year  of  statistically  prepared  STATS 
data  to  produce  a modified  statistical  file  for  use  with  the  DSPA  program. 
The  procedure  involves  the  scaling  of  the  STATS  file  data  on  a monthly 
basis  by  factors  computed  from  user-specified  proportions  and  confidence 
levels.  In  addition,  the  PROFILE  program  will,  upon  request,  perform  a 
worst  case  analysis  for  low  solar  insolation,  low  wind,  or  high  wind 
periods.  The  exact  type  and  use  of  these  input  requests  and  their  results 


are  discussed  below.  The  revised  weather  data  is  written  onto  a STAT 
file  for  subsequent  use  as  input  to  the  DSPA  program. 

1 ) Low  or  High  Temperature  Months 

In  constructing  the  hour-by-hour  profile  year  for  temperature,  the 
following  options  are  available  for  each  month; 

a)  Use  the  average  hourly  values  for  temperature  direct.,  from  the 
STATS  file. 

b)  Construct  a month  for  which  the  temperatures  are  low,  in  the 
sense  described  below,  at  a user-specified  confidence  level. 

c)  Construct  a month  for  which  the  temperatures  are  high,  in  the 
sense  described  below,  at  a user-specified  confidence  level. 


A low  temperature  month  is  defined  as  one  having  an  average  tempera- 

low  Hi 

ture  C„  such  that  100P%  of  all  actual  monthly  average 

Urn,  temp.  ! . 

temperatures  can  be  expected  to  exceed  temp  ^ lOOU/o  confi- 

dence level.  The  proportion  P and  the  confidence  level  a are  input 

low 

by  the  user.  The  PROFILE  program  calculates  C„  ^ using  P,  a, 

im,  temp. 

the  monthly  temperature  average  for  each  year  of  data  (statistic  Ua 

of  Section  3.3)  and  the  mean  and  standard  deviation  of  these  values 

(5a  and  6a  of  Section  3.3).  The  difference,  . = f.  - 

’ Jim,  temp.  Jim 

C„  is  used  to  scale  down  all  hourly  temperature  values  for 

Jim,  temp. , 

the  month. 


Similarly  a high  temperature  month  is  defined  as  one  having  an 

average  temperature  such  that  lOOP^  of  all  actual  monthly 

Jim, temp.  ... 


temperatures  can  be  expected  to  fall  below  C 


confidence,  and  the  difference. 


.high 
im,temp, 


phigh  _ 

,temp  Jim’ 


- C 


with  100a$ 
high 


Jim, temp. 

used  to  increase  each  hourly  temperature  value  for  the  month. 


IS 


2)  Low  or  High  Solar  Insolation  Month 

The  options  and  descriptions  of  the  low  or  high  solar  insolation 
construction  are  exactly  analogous  to  those  for  temperature  above, 
except  that,  since  solar  insolation  is  on  a ratio  scale  (i.e.,  has 
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a naturally  occurring  zero  point  and  equal  interval  measurements),  a 


multiplicative  scaling  factor  is  more  appropriate  than  the  differ- 

low 

ence  technique  applied  to  temperature  scaling.  Hence, 


,low 


,high 


_ phigh 


iva,  insol. 


„ . ana  n„  ^ . , / Q„  , where  Q.  is  the 

£m,  insol./  Jim  Jim,  insol.  Jim,  insol./  Jim  Jim 


C„:  ' / Q„  and  R 

observed  monthly  solar  insolation  (statistic  5c  in  Section  3.3). 


3)  Low  or  High  VJind  Velocity  Month 


The  explanations  and  computations  involved  in  scaling  the  hourly 

values  of  wind  velocity  to  produce  a low  or  high  month  are  the  same 

low 

as  those  for  solar  insolation.  The  scaling  ratios  R.  • j i 
...  Jim,  wind  vel. 

and  R„  ® , are  used  to  lower  or  raise,  respectively,  the 

Jim,  wind  vel.  > r 

hourly  wind  velocity  observations. 


U ) Period  of  Low  Solar  Insolation 


In  addition  to  producing  uniformly  low  or  high  or  mean  solar 
insolation  values  for  any  given  month,  the  PROFILE  program  also 
provides  an  option  for  creating  a string  of  "worst  case"  days  for 
which  the  daily  insolation  values  fall  below  some  user-specified 
fraction  of  the  observed  average  for  that  day.  The  fraction  P and 
probability  of  occurrence  a are  input  by  the  user  for  each  month. 

For  each  month,  in  which  a low  period  is  requested,  the  PROFILE 
program  compares  the  daily  average  for  a given  day  of  each  data  year 
(statistic  2a  of  Section  3.3)  with  the  average  over  all  years 
(statistic  3a).  If  the  average  for  a particular  day  is  less  than 
lOOP^  of  the  average  for  the  same  day  over  all  remaining  years, 
then  the  day  is  classified  as  "low";  otherwise,  :he  day  is  termed 
"normal . " 

A pattern  of  "low"  and  "normal"  days  is  thus  created  for  a given 
month  for  all  data  years.  The  lengths  an<i  the  number  of 

occurrences  (n)  of  strings  of  "low"  days  are  considered  as  realiza- 
tions of  a random  variable  having  a Poisson  type  distribution  (rare 
event).  The  PROFILE  program  uses  the  patterns  of  these  "low"  strings 
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to  determine  the  average  number  of  consecutive  "worst  case"  days 

The  value  of  as  discussed  in  Appendix  B,  is  a Justifiable 


Jim 


approximation  to  the  true  Poisson  maximum  likelihood  estimate, 

X.  From  the  value  of  (hereafter  referred  to  as  X)  and  the  user- 

x-m 

specified  confidence  level  a,  the  program  computes  a string  length 
N such  that  the  probability  of  N or  fewer  "low"  days  in  a rev; 
is  a. 


This  N-day  string  is  centered  about  the  15th  day  of  the  month,  aid 
the  hourly  solar  insolation  values  for  each  of  these  days  is  scaled 
down  by  the  factor  P.  It  is,  of  course,  possible  that  a user- 
defined  "low"  day  may  be  higher  than  the  same  day's  values  in  the  low 
solar  insolation  month  generated  by  item  2 above.  In  such  cases, 
the  PROFILE  program  will  select  the  minimum  of  the  two  computed 
values. 

The  number  of  sequential  "low"  days  generated  is  directly  determined 
by  the  values  of  a and  X.  To  assist  the  user  in  planning  "worst 
case"  insolation  analyses.  Appendix  D contains  a table  of  the 
number  of  "low"  days  as  a function  of  a and  X.  The  value  of  X is 
dependent  on  the  particular  station,  month  of  the  year,  and  user- 
input  fraction  P.  Appendix  C gives,  for  each  station  (see  Table  1-1 
for  a cross  reference  to  station  location  and  station  number),  a 
table  of  X (actually  as  a function  of  the  month  and  P.  Hence, 

the  user  may  select  a P and  an  a for  a particular  station  and  month 
and  determine  the  resultant  PROFILE  program  X from  Appendix  C,  and 
subsequently,  the  resultant  number  of  sequential  "low"  days  from 
Appendix  D.  The  user  may  then  vary  either  P (to  get  a different 
X)  or  a to  obtain  a more  desirable  "low"  string  length. 

Note  that  the  largest  value  of  X found  is  8.235  for  P = 1.00  for 
month  7 of  location  93193.  Using  this  value  for  illustrative 
purposes.  Appendix  D indicates  that  a string  of  13  "low"  solar 
insolation  days  will  occur  for  a = 0.93  (93^).  However,  a P of  1.00 
is  really  somewhat  useless  for  a "worst  case"  since  all  values 
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would  remain  as  100^  of  the  original  averages.  More  reasonably,  for 
P < 1.00,  the  value  of  X will  be  less  than  2.0  yielding  "low" 
strings  of  fewer  than  6 days  for  Qi  < 0.99  (99%).  In  any  case,  with 
the  aid  of  Appendix  C and  Appendix  D,  the  user  can  predetermine  the 
results  of  his  "worst  case"  analysis. 

5 ) Period  of  Low  or  High  Wind  Velocity 

As  in  the  case  of  low  solar  insolation  periods  (item  1*  above),  a low 
or  high  wind  velocity  day  is  defined  as  one  in  which  the  average 
daily  wind  velocity  is  less  than  lOOP^^  or  greater  than  lOOP^^  of 
the  dally  average  for  the  same  day  over  all  years.  The  PROFILE  pro- 
gram fits  a Poisson  distribution  to  the  strings  of  low  days  and  of 
high  days.  The  user-input  occurrence  probabilities,  a^,  and  are 
then  used  to  determine  the  number  of  sequential  "low"  days,  N^,  and 
the  number  of  sequential  "high"  days,  Ng,  in  the  same  manner 
described  for  low  insolation  periods  above. 

The  "low"  days  are  centered  around  the  10th  day  of  the  month,  and 
the  Ng  high  days  are  centered  about  the  20th  day.  The  "low"  days 
are  then  scaled  down  by  the  factor  P^,  and  the  "high"  days  are  scaled 
up  by  the  factor  (>  l.OO).  The  remaining  days  of  the  month  may 
be  taken  from  the  mean  year  data  or  from  low  or  high  month  values. 

The  mathematical  formulas  used  for  all  of  the  above  PROFILE  computations 
are  provided  in  Appendix  B.  A flow  diagram  of  the  PROFILE  program  is 
provided  in  Figure  3-5. 
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FIGURE  3-5.  PROFILE  COMPUTER  PROGRAM 
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APPENDIX  A 

STATS  PROGRAM  MATHEMATICAL  FORMULAS 
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Input  (frcm  Merge  Tape): 


'^Imdhy  ’ 


, ,,  j Q-,  = Temperature,  wind  velocity,  solar  insola- 

^ my  tion  for  location  1,  month  m,  day  d,  hour  h, 
year  y 


Y,  = Humber  of  years  of  data  for  month  m,  at 
Im  . 

location  X 


D = Number  of  days  in  month  m 
m 

= Total  number  of  years  of  data  at  location  1 


Statistic  Formula 


la) 

T 

Imdh 

= 

Y't  /y 

^ Imdhy  Im 

lb) 

V 

Imdhy 

= 

y V /Y 

Y Imdhy Im 

; 

( 

Ic) 

*^lmdh 

= 

^ '^Imdhy'^^lm 

2a) 

V 

Imdy 

= 

^ ^Imdhy^^^ 

2b) 

^Imdy 

= 

^‘^Imdhy'^^^ 

■ 

3a) 

"imd 

= 

EVdy/^lm 

3b) 

“^Imd 

= 

^ ‘^Imdy'^^lm 

Ua) 

T-, 

Imy 

= 

/2UD 

d h 

V 

i 

Ub) 

V, 

Imy 

= 

23  23  V,  ,,  /2bD 
^ “ Imdhy  m 

d n 

he ) 

Q 

^Imy 

= 

23  23  Qt  /2bD 

" ” imdhy  m 

d h 
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PROFILE  PROGRAM  MATHEMATICAL  FORMULAS 
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1.  Tolerance  Limits  - Temperature. 


Input;  a,  P,  T^^,  ST^^ 


za  = (1  - a),  z^_p  = (P) 


where:  $(x)  = 


^ -t^/2 

J 


a,  = 1 - z^  / 2(Y,  - 1) 

Im  a / |_  Im 


hm  “ 4-F  - 4'hm 


Im 


=1-P  " V'Lp  - 


a b 
Im  Im 


Im 


aow  = T,  - k,  ST, 

Im,  temp.  Im  Im  Im 


high  ^ .j,  ^ g.p 

Im,  temp.  Im  Im  Im 


2.  Tolerance  limits  - solar  insolation. 


. = Q _ k,  SQ- 

Im,  insol.  Im  Im  Im 


-high  — Q + k *^0 

^Im,  insol.  " ^Im  *^lm  “^Im 
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3.  Tolerance  limits  - wind  velocity. 


^low 


c:“"  . , , = V,  - k,  SV, 

Im,  wind  vel.  Im  Im  Im 


chigh  = V.  + k,  SV, 

Im,  wind  vel.  Im  Im  Im 


1*.  Translation  Factors  for  low/high  temperature  months. 


plow  ^ - _ plow 

Im,  temp.  Im  Im,  temp. 


phigh  _ ^ plow 

Im,  temp.  Im  Im,  temp. 


5.  Scaling  Factors  for  low/high  solar  insolation  months. 


plow  ^ plow  / ^ 

Im,  insol.  Im,  insol./  Im 


high  ^ / Q, 

Im,  insol.  insol.  J Im 


6.  Scaling  Factors  for  low/high  wind  velocity  months. 


plow  ^ plow  /y 

Im,  wind  vel.  Im,  wind  vel./  Im 


phigh  ^ (.high  hj 

Im,  wind  vel.  Im,  wind  vel./  Im 
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Generation  of  period  of  low  solar  insolation  days. 


Input:  P,  a,  Q,  , Q,  , 

Imay,  Imd 


Im 


1 if  Q.  , < P • Q?  ^ 
tody  — Jcmd 


“ “■  “urndy  » P • “U 


where  is  the  mean  of  all  values  exlcudin^  the  current  one. 

A value  of  X,  is  recorded  for  each  day  of  data  for  month  m.  Patterns  of  O's 
im 

and  I's  are  thus  observed. 


Let  = n if  the  pattern  0,  1,  1,  . . . , 1,  0 is  observed,  where  there 

are  n consecutive  I's.  The  mean  B,  of  all  B,  values  is  calculated  and  used 

Im  Im 

— 

as  X (in  Reference  8)  to  produce  an  estimate  X of  the  parameter  X in  a truncated 
Poisson  distribution. 


The  truncated  Poisson  distribution  is  summed  until  the  cumulative  proba- 
bility level  is  greater  than  a,  that  is,  until  a value  li  is  produced  such  that 


N 

x:xVi! 

D 

m .. 

Z^Vi! 


< a 


and 


W+1 


x:xVi! 

D 

m 


> a 
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ITie  N low  solar  insolation  days  are  produced  having  hourly  insolation 
values : 


0 = P • 0 

^Imdh  ^Imdh 


8.  Generation  of  periods  of  low/high  wind  velocity  computations  are  exactly 
the  same  as  in  except  that  2 strings  are  generated,  one  for  low  wind  velo- 
city days  (P^  < l)  and  one  for  high  wind  velocity  days  (P^  > l).  low  days 
having  hourly  wind  velocities; 


low 

Imdh 


V 


Imdh 


and  Ng  high  days  having  hourly  wind  velocities : 

^high  _ p . — 

Imdh  2 Imdh 

are  then  generated. 

Distribution  Fitting  in  the  PROFILE  Program 

The  Profile  program  calculates  the  average  number  of  consecutive 

low  days  observed  for  location  i,  month  m in  all  years. 

The  distribution  to  be  fitted  is  a right-truncated  Poisson  with  parameter 

X.  The  truncation  point  is  D , the  number  of  days  in  month  m.  D = 28,  30  or 

m m 

. 31.  The  frequency  function  is 

! 

S 

[ f(x;  X,  D„)  - KX-Zxl  X = 0 


where 
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In  the  PROFILE  program,  we  need  to  calculate  the  (l  - a)  tail  of  this  distri- 
bution. We  need  an  integer  D*  such  that 


i.i 

y ] f(xi  X,  D^)  s 1 - a 


x=D*» 


To  do  this,  an  estimate  of  the  unknown  parameter  A is  needed.  The  usual 


technique  of  maximum  likelihood  produces  an  estimate  X which  is  the  unique 


.th 


positive  real  solution  to  the  degree  polynomial  equation 


/ TT 

V'  I ^ 


X^  = 0 . 


(1) 


Reference  8 shows  that  when  itself,  rather  than  the  maximum  likelihood 
estimate  X,  is  used  to  estimate  X,  the  error  is  negligible  for  values  of 
and  encountered  in  the  PROFILE  program.  For  derivation  of  equation  (l) 
and  more  details  on  accuracy  of  B£^  as  an  estimate,  see  Reference  8. 
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